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ABSTRACT 

In this paper we examine the influence of stochastic effects on the integrated 
near-infrared hght of star clusters with ages between 7.5 < logt < 9.25. To do 
this, we use stellar evolution models and a Monte Carlo technique to simulate 
the effects of stochastic variations in the numbers of main sequence, giant, 
and supergiant stars for single-generation stellar populations (SSPs). The 
fluctuations in the integrated light produced by such variations are evaluated 
for the VJHK bands. We show that the V JHK light of the star clusters can 
be strongly affected by plausible stochastic fluctuations in the numbers of bright 
but scarce stars. In particular, the inclusion of thermally pulsing AGB stars in 
the stellar evolution models yields integrated colors with values in agreement 
with the spread seen for Large Magellanic Cloud clusters that are known to have 
signiflcant number of AGB stars. Implications of this analysis are important for 
studies of the integrated light of stellar populations where it is not possible to 
resolve individual stars. 



Subject headings: galaxies: star clusters — galaxies: stellar content — infrared: 
stars — methods: statistical — stars: evolution 



1. Introduction 

The main question addressed in this work is how large is the contribution of small 
numbers of bright stars to the near-infrared (near-IR) integrated light of a star cluster? 
For a cluster of given age and metallicity, the size of the contribution will depend on 
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the evolutionary state of the stars in question as characterized by their temperature and 
luminosity. More generally, the effects of changing small numbers of stars will depend on 
the total mass of the cluster and the cluster's initial mass function, its IMF. For stochastic 
variations of a given size the integrated light from less massive clusters will obviously be 
more strongly influenced than will the light from more massive clusters. The advantages of 
including the JHK bands in the analysis is that the effects of variations in the numbers 
of later type giants and supergiants in a cluster on the integrated light will be much more 
evident than if only optical wavelengths were examined. 

The age range analyzed (7.5 < logt < 9.25) was chosen because the available isochrones 
contain all relevant stellar evolutionary phases up to the end of the C-burning phase for 
stars with M > TM© and up to the end of the early asymptotic giant branch (E-AGB) 
for stars with 2 < M/M© < 5; for stars with M < 1.7M0, corresponding to logt > 9.25, 
the evolution is followed up to the He-flash (pchaller et al. 1992 ). Consequently, although 



the thermally pulsing asymptotic giant branch (TP-AGB) for intermediate mass models 
is missing, the He-burning phases are covered by the models in the age range considered. 
The TP-AGB according to models by [Bedijn (1988) was taken into account in a later stage 



of this work providing information on its contribution to the integrated light of stellar 
populations and allowing a better approximation to real star clusters concerning their 
integrated colors and stochastic fluctuations. The lower age limit selected avoids massive 
stars in evolutionary stages with very high mass loss rates, for which bolometric correction 
and temperatures are very uncertain. The present work is restricted to solar metallicity 
since the bolometric corrections and the color-temperature relations available come from 
solar neighborhood stars. This restriction, though, do not affect our primary goal which 
is to examine stochastic effects in stellar populations on the integrated colors of those 
populations. 



Similar studies using visible bands have been carried out by, e. g., [Barbaro fc BerteHi 



(1977), 


Becker k Mathews (1983) 


, Uhiosi et al. (1988) 


al. (1995) 


. In particular, |Chiosi et al. (1988) employ a 



models that include convective overshooting and cover the whole range of evolutionary 
phases. They calculated integrated {U — B) and [B — V) colors and synthetic H-R diagrams 
of star clusters as a function of age, and applied their results to the known age and color 
distribution of clusters in the Large Magellanic Cloud (LMC). Beside the stochastic effects 
on the IMF, they also take into account the effects of dispersion in the stellar ages. In 
summary, among other important results, their analysis indicates that integrated UBV 
colors show a significant spread for ages at which few luminous AGB stars are present. 
Thus, they conclude that it is necessary to take stochastic effects into account when the 
integrated light of a cluster is used to estimate the cluster's age and chemical composition. 
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pirardi et al. (1995)| employ the photometric models of simple stellar populations (SSPs) 
by [Bertelli et al. (1994)| and the integrated UBV colors of LMC star clusters ( [Bica et "aT] 
|1996D to investigate the history of star formation in the LMC. They present theoretical 
integrated UBV colors of SSPs with ranges in age and metallicity that span the observed 
ranges for the clusters. They analyze the dispersion in the observed integrated colors of 
the clusters and conclude that in UBV stochastic effects dominate the color dispersion of 
the young clusters. Dispersions in the colors due to variations in reddening and metallicity 
become significant for clusters older than 300 Myr. In summary, the highest dispersion in 
the integrated colors of the SSPs occurs for ages between 6.3 and 25 Myr as a consequence 
of the number and luminosity of red supergiants. Similar results were obtained for different 
metallicities, except for smaller dispersions in young, metal poor SSPs where red supergiants 
are rare compared to a more metal rich population. 

The near-IR light of a stellar system can be dominated by only a small number of 



luminous, red stars most of whose energy is emitted in the near-IR (e. g. Persson et al 



19831 , [Frogel et al. 19901) . Therefore, stochastical variations in the numbers of these stars 



could lead to erroneous conclusions concerning the age and/or metallicity of the system if 
such estimates are based on integrated light measurements. Thus, in the present study we 
analyze the effect on the integrated JHK light from a single-generation stellar population, 
of stochastic variations in its IMF. This issue is addressed here in the hope of setting useful 
constraints on stellar populations of galaxies, allowing one to decide how much an integrated 
observation in JHK of a star system of certain age and metallicity can be affected by the 
stochastic nature of the IMF in single-generation stellar populations. 

In section 2 the data employed and the procedure carried out to transform the 
theoretical H-R diagram into observational color-magnitude diagrams are described. 
Section 3 presents the Monte Carlo technique used to carry out the simulations and a 
comparison between a young SSP and an intermediate age one. Section 4 discusses the 
relationship between the stochastic effects and the age of a SSP. In section 5 the results 
are compared with observations of LMC star clusters. Finally, section 6 presents the 
conclusions. 



2. Data 
2.1. Stellar Evolution Models 



The stellar evolution models by Schaller et al. (1992, hereafter SSMM) were employed 
in the present study. Basically they present models from 0.8 to 120 M©, at metallicities 
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Z = 0.001 and solar, with mass loss and convective overshooting. The stellar tracks were 
transformed to isochrones by using the code made accessible by G. Meynet. The fainter 
end of the main sequence (MS) (m < O.SM©) was extended down to O.OTM© by using the 
information in [Bessell (1991)| (Teg and L) and [VandenBerg et al. (1983)] (mass-luminosity 
relation to get m). For most of the ages investigated less massive stars have not yet 
contracted to the MS phase, as at least 1 Gyr is needed for a complete MS to be formed 
( [Mera et al. 1996] , and references therein). However, a complete MS was used here, even for 



the younger SSPs, as a rough approximation to real clusters where low mass stars are in 
pre-MS phases. This approximation should be a lower limit in terms of near-IR integrated 
light contribution, since low mass contracting stars are cooler and brighter than their future 
MS location. 

In the present context, the fact that the models do not take into account evolution 



beyond the E-AGB for intermediate mass stars must be considered. [Ferraro et al. (1995) 
carried out BVJHK observations of stars in 12 intermediate age LMC clusters showing 
that around 60% of their integrated K light comes from AGE stars, consistent with the 



results of Persson et al. (1983) and Frogel et al. (1990). [Chariot fc Bruzual (1991 



constructed stellar population models by using [Maeder fc Meynet (1989)| stellar evolution 



models complemented with observations of TP-AGB. Their models indicate that the whole 
AGB contribute with 60% at J and 70% at K to the integrated light of a 0.4-2 Gyr SSP, 
the TP-AGB stars alone contributing with 20% at K and a little less at J. The agreement 
between these studies and their results support the conclusion that the E-AGB is the major 
contributor to the whole AGB in the integrated near-IR light of intermediate age star 
clusters. The TP-AGB phase, though, has not negligible contribution and should be kept in 
mind, especially when stochastic variations in such stellar populations are to be evaluated. 



2.2. Near-Infrared Colors 



The JHK (and V , in the Kron-Cousins system) colors, in a standardized system, are 
from [Bessell fc Brett (1988)| for MS stars from B8 to K5 and for giants from GO to M6, and 
from Bessell (1991, hereafter B91) for later MS spectral types. When needed, the colors 
for MS stars earlier than B8 were linearly extrapolated using the hottest spectral types 
available; the same has been done for giants earlier than GO. Colors for 05-M5 supergiants, 
all considered as evolved stars with Mhni < —4, are from [Koornneef (1983)[. In this case, 



the transformations in [Bessell fc Brett (1988)[ were used to bring Koornneef 's colors into 
their system. 
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2.3. Bolometric Corrections and Temperature Scale 

In order to transform the theoretical H-R diagram into a color-magnitude diagram 
(CMD) for comparison with observations, bolometric corrections and temperature scale 
compiled by Schmidt-Kaler (1982, hereafter SK82) for giants, supergiants and for MS stars 
of spectral type earlier than M and those determined by B91 for fainter stars were employed. 
The link between the colors and BCy and Tefj was guided by the spectral type provided in 
both sets of data. An exception to this procedure was adopted for the M dwarfs, for which 
B91 gives colors, Tes, and all other information needed to get BCy (his Table 2). 

The bolometric correction was defined here as BCy = Mboi — My and its zero-point 
was set by adopting BCy© = —0.08 (Mboi© = 4.75). The consistency between B91 and 
SK82 data was checked by looking at the flux density of a zero-magnitude star in both 
scales. The V magnitude difference at ^ = is 0.03 mag brighter in the Johnson system 
(adopted by SK82), which is insignificant given the uncertainties. The temperature scale 
adopted was checked at spectral type K7V (where the data from the different sources were 
merged) and no significant discrepancy was found. On the other hand, at this spectral type 
there was a difference of 0.09 mag in the bolometric correction of SK82 and B91 (in the 
sense that BC is larger in SK82). The option adopted was to smooth out the BC around 
Tefj = 4000 K (i.e. K7V) as long as this procedure is confined to the fainter end of the MS 
and will not have a major effect on integrated light models. 



3. Simulations of Single-Generation Stellar Populations 



[ 



] 



Barbaro & BerteUi (1977), [Becker fc Mathews (1983)^ [Chiosi et al. (imj , |Girardi 
fc Bica (1993)1 and [Girardi et al. (1995)| investigated the importance of stochastic effects 
on the IMF of star clusters in the UBV bands. Similarly, in the present study we used a 
Monte Carlo technique as an approach to assess the impact of stochastic variations in the 
IMF of single-generation systems on near-lR bands. 

The Monte Carlo integration used simulates the stellar mass distribution of a system of 
a given age: a randomly generated mass is weighted on the IMF and the process is repeated 
until the number of stars sampled reachs a preset value. A total of 100 simulations of 
systems containing a variable number of stars (10^ to 10®) is performed and the mean and 
dispersion of the integrated colors are computed. 

Following is a more detailed description of the procedure used to carry out the 
simulations. An introduction to Monte Carlo methods is given by, e. g., [Allen et al. (1996] 
and Press et al. ( 1992 
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3.1. Theoretical Background 

The IMF of a steUar system gives the relative number of stars within a mass range and 
it is employed in its usual definition 

$(m) = dN/dm = Am-(^+^\ (1) 

where A is a constant and x the IMF slope, x = 1.35 is the value corresponding to 
the [Salpeter (1955)| luminosity function. A is determined by integrating the IMF from the 



lowest (m;=0.07) to the highest stellar mass (mu=120) in the system and by normalizing 
the number of stars in the system to 1. Specifically, 

^ = ^ (2) 

This normalization is critical for the present application of the Monte Carlo method since 
now $(m) can be treated as a probability distribution function (PDF) which gives the 
probability that a random mass m' is in the range between m and m + dm. $(m) now 
obeys the following conditions: 

$(m) > 0, for any m (3) 

u 

$(m')rfm' = 1 (4) 

'nil 

where mi and m„ are the lower and upper mass limits of the stellar system. 

Now, the PDF $(m) can be transformed into a new PDF g{N), where is a 
single-valued function of m. Therefore, the probability of occurrence of the random variable 
N' within dN and the probability of occurrence of the random variable m' within dm must 
be the same, i. e., 

|$(m)dm| = \g{N)dN\ (5) 

From Eq. (p, a cumulative distribution function (CDF) can be constructed which 
gives the probability that the mass m' is less than or equal to m, i. e., 

N{m) = / ^m')dm' (6) 

J nil 

and, consequently, dN/dm = $(m), resulting from Eq. (|^) 

g{N) = 1, < iV < 1 (7) 



This means that the CDF is a uniform function for which any value is equally likely in 
the interval [0,1]. 
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The CDF property of uniformity between and 1 allows one to use a random number 
generator (namely, that of Parker and Miller with Bays-Durham shuffle as presented by 
Press et al. (1992}] ) to sample values for the CDF N{m). Since the random variable m and 
the CDF N{m) are related uniquely, it is possible to sample m by sampling first and 
then to solve the result of the integration of Eq. @ for m, i. e., 

m = [mY%l -N)+ m-^'AT]-!/^ (8) 

Each stellar mass m generated in this way is interpolated linearly on an isochrone in order 
to get the corresponding log Teg and logL. The correlations assembled in section 2.3 were 
then linearly interpolated to the colors corresponding to these values of logTgfj and logL. 
The mass difference between consecutive points on an isochrone is proportional to the 
evolutionary lifetime. Therefore a random mass will have larger probability to be sampled 
in evolutionary stages for which there is a relatively wider range of masses indicating longer 
lifetimes. Consequently, for a SSP of a given age, the fact that the number of stars in each 
evolutionary phase is proportional to its lifetime is taken into account. Integrated fluxes 
and colors are then computed by appropriately combining the fluxes and colors from the 
individual stars with the weights as described above. 

3.2. A Comparison Between Young and Intermediate Age SSPs 

In this section, a comparison of the simulations carried out for solar metallicity SSPs 
with ages 32Myr (logt = 7.5) and 1 Gyr (logt = 9.0) is presented. With the adopted 
criterion that stars with Mboi < —4 are supergiants, it turns out that the post MS stars in 
the 32Myr SSP are slightly brighter than this limit. We also computed simulations treating 
these stars as giants for comparison. Although the colors of the giants and supergiants 
differ, we flnd no signiflcant differences between these two sets of simulations apart from 
those resulting from the stochastic effects. 

Fig. 1 shows how the young stellar population is distributed in CMDs formed by both 
visible and near-IR bands (a,c) and exclusively by near-IR bands (b,d). Figs, la and lb 
show different CMDs for the same SSP simulated, while Figs. Ic and Id show the same 
for another simulation. Each SSP has 10^ stars and an IMF with x = 1.35. Only the 
upper MS and the supergiants, the only regions for which the Monte Carlo simulations 
reveal signiflcant variations, are shown. The integrated colors corresponding to each CMD 
are indicated on the panels. Fig. 2 corresponds to Fig. 1 for an intermediate age stellar 
population. We also ran an experiment to see, in the absence of any stochastic effects, how 
many stars would be required to produce at least one post MS star in such SSPs. We found 
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that young SSPs (logi = 7.5) containing 140 stars for an IMF slope x — —0.5, 8990 stars 
for X = 1.35 and 136540 stars for x — 2.0 should contain one post MS star if a continuous 
IMF is considered. Similarly, intermediate age SSPs (logt = 9.0) containing 170 stars for 
an IMF slope x = —0.5, 890 stars for x = 1.35 and 5590 stars for x = 2.0 should contain 
one post MS star. For SSPs less populous than these the stochastic nature of the IMF can 
still produce post MS stars and, consequently, color variations as large as those found for 
SSPs with number of post MS Npms=l. If an SSP is far less populous then there would be 
no post MS stars, and the integrated fluxes would be dominated by upper MS stars with 
little spread in color, which results in smaller color fluctuations. 

Fig. 3 presents the integrated colors (V — K), {J — K), {H — K) as a function of 
the number of post MS stars for the young SSP (a, b and c) and the older one (d, e and 
f), where 3 slopes of the IMF are also considered. Fluctuations in the colors come from 
100 simulations for SSPs of fixed total number of stars and are represented by twice the 
standard deviation from the mean color. {V — K) has the largest fluctuations because it is 
most sensitive to the slope of the spectral energy distribution of the stars in the system. 
The fluctuations are also larger the fewer post MS stars there are in the system. As their 
number increases, not only do the color fluctuations decrease but the mean colors get redder 
as well. As may be seen from Fig. 3, the stochastic effects are always more signiflcant than 
variations due to changes in the IMF slope for the range studied here except for the young 
SSP when the number of post MS stars gets to be 100 - then {V — K) is larger for x — —0.5 
than it is for steeper values of the IMF slope. 

Relative flux contributions from MS and post MS stars for well-populated SSPs with 
100 post MS stars are plotted for the 3 slopes of the IMF in Fig. 4. For the young SSP 
(left side of the Fig. 4) the V flux contribution is distributed between the MS turnoff (B 
type stars) and the supergiantes, with insignificant contribution coming from the lower 
MS. As X increases from -0.5 to 2, the MS turnoff V contribution increases from 40% to 
60%, while the supcrgiants one decreases from 60% to 30%. The supergiants dominate the 
integrated JHK bands for all IMF slopes, followed by the MS turnoff, which contributes, 
at most, with 15% of the near-IR light. For the older SSP (right side of Fig. 4) the V flux 
contribution is distributed again between the MS turnoff (types A and F) and the post 
MS stars. As x increases from -0.5 to 2, the turnoff (A+F stars) V contribution increases 
from 53% to 65%, while the giants one decreases from 45% to 30%. Differently from the 
young SSP, there is a significant contribution in the near-IR bands from the turnoff: as x 
increases from -0.5 to 2, A and F stars together increase their contribution to H and K 
from 20% to 27%, the contribution in J reaching 40%. The giants contribute with 80% 
(x = -0.5) to 60% {x = 2) in H and K and with 72% {x = -0.5) to 53% {x = 2) in J. 
We can compare these results with that found for an old (~ lOGyr) stellar population (e.g. 
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brogel 19"88D . For such a population, evolving dwarfs, subgiants and K giants dominate the 



B and V integrated light, regardless of the IMF. In JHK the M giants are the largest single 
contributors to the flux. For x = 1.35 their contribution reaches 50% of the total K band 
light and is greater still for x = 0. In connection with the present study, Frogel's models 
show the relative importance of red giants in the integrated light of composite systems that 
contain an older stellar component, and whose contribution can change from 1/ to by 
significant amount. 

Subsequent analysis (section 5) shows that the inclusion of the TP-AGB can influence 
significantly the integrated colors in the case of the SSP with logt = 9. 



4. SSP Age and the Stochastic Effects 

For the whole range of ages (7.5 < logt < 9.25) SSPs were computed in steps of 
Alogt = 0.25. The behavior of the integrated (J — K) color as a function of logt for SSPs 
with different total masses and IMF slope x = 1.35 is shown in Fig. 5. Panel (a) shows 
the color evolution of SSPs with 10^ and 10^ stars; it is also shown the percentage of post 
MS stars for each age (=Npms/(N.100)). When SSPs of identical total number of stars are 
compared, they can be associated to the evolution of a single star cluster. However, young 
clusters with a small number of stars are more difficult to interpret in terms of evolution 
since their post MS may be underrepresented, which produces large fluctuations in the 
integrated colors. Fig. 5a shows that: 

• Less massive clusters with 10^ stars present a gradual reddening as they age. 

• Young massive clusters with 10^ stars can have integrated (J-K) as red as those of 
older clusters. 

• Young star clusters have integrated near-IR colors very sensitive to their masses. 

• Older star clusters of different masses have small spread in color. 

In Fig. 5b the color evolution is presented for SSPs with a fixed number of post MS 
stars. The main difference between systems with Npms=100 and Npms=500 is the smaller 
color fluctuations for the last ones, the mean color having similar behavior to that of the 
more massive SSPs with N= 10^ stars (see Fig. 5a). 

Well-populated SSPs with at least 100 post MS stars have their integrated colors 
modulate by a combination of factors like the ratio between the number of upper MS 
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and post MS stars and the extension in temperature of the blue loop and its luminosity 
compared to the asymptotic giant branch and to the MS turnoff. 

For X = 1.35, color fluctuations of 10% (la from the mean) are reached when there are 
70 post MS stars for logt = 7.5, 80 for logt = 8.0, 350 for logt = 8.5 and 40 for logt = 9.0. 
The larger the number of post MS stars, the smaller the color fluctuations. The extended 
E-AGB at age logt = 8.5 is the main feature responsible for the largest fluctuations in the 
integrated color of that SSP compared to the other ages. 

Fig. 6 presents the evolution of the VJHK flux contribution of MS and post MS stars 
to the integrated light of SSPs built with x = 1.35 and 100 post MS stars. In summary, the 
results indicated: 

• In V, as the age increases from logt = 7.5 to 8.25 the MS B stars contribution gets 
between 55% and 63%, while the post MS contribution decreases from 40% to 25%. 
The contribution of redder MS spectral types being less than 15%. 

• For logt = 8.5 MS B, MS A stars and post MS stars share most of the total integrated 
V light, contributing with 40%, 27% and 23% respectively. The major contributors 
for the V light of older SSPs change dramatically: for 8.5 < logt < 9.25, the post 
MS contribution increases from 25% to 35%, the MS A stars decrease strongly their 
contribution from 60% to 0% while the MS F ones increase theirs from 10% to 55%. 

• In H and K, the post MS contribution dominates throughly the near-IR integrated 
light from ^90% (at logt = 7.5) to ^70% (at logt = 9.25); the contribution of the 
turnoff increases for the older SSPs reaching 25% at most. 

• In J, the turnoff MS types have signiflcant contribution especially for the older SSPs: 
the B type peaks at logt = 8.25 with 23%, the A type at logt = 8.75 with 32% and 
the F type at logt = 9.25 with 30%. 



5. Comparison with Observations 

We compare the integrated color models and their fluctuations to the integrated colors 
of Large Magellanic Cloud clusters (from Persson et al. 1983| ) in Fig. 7a. Regarding the 



age spread of the models computed, they are in principle adequate for a comparison with 
clusters of SWB types ranging from II to V. However, in most of the cases, the models are 
not expected to encompass clusters with SWB type IV and V because the stellar evolution 
models do not take into account evolution beyond the E-AGB for intermediate mass stars. 
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where bright AGBs hke Carbon stars can contribute significantly to the integrated hght of 
a cluster. 

Frogel et al. (1990) analyzed the AGB contribution to the integrated light of Magellanic 
Cloud clusters. They show that Carbon stars (originated by TP-AGB stars) are only found 
in clusters of SWB types IV- VI and account for 50% to 100% of the bolometric luminosity 
from the AGB above the RGB tip (He-fiash). The same type of stars contribute with 30% 
to the bolometric luminosity of a cluster with age around 2 Gyr (SWB V-VI), the whole 
AGB (Mboi < -3.6) contributing with 40%. 

The agreement between the younger SSP models with clusters of types SWB II and 
III in the color-color diagram of Fig. 7a is indicating that the isochrones are essentially 
complete for ages around logt = 7.5. More importantly, the small scatter of clusters 
around the younger SSP models indicates that their integrated colors are compatible with 
SSPs showing an almost continuous IMF, since big fiuctuations in color are not needed to 
explain the young cluster distribution. It seems that the inclusion of the TP-AGB phase in 
the stellar evolution models would improve the older SSP integrated colors, matching the 
observations. 

In order to address the aforementioned point, the infiuence of the TP-AGB on the 
integrated colors was evaluated quantitatively by complementing the isochrones with the 
upper AGB semi-empirical models of Pedijn (1988)| , which consider that stars with initial 
masses in the range 0.8 < M/M© < 8 will pass through that phase. The BCj^ vs (J — K) 
and Mboi vs log T^s relationships in [Frogel et al. (1980]| together with the AGB tip near-IR 



colors for LMC clusters observed by [t'erraro et al. (1995)| were used. The TP-AGB has 



a maximum extension for the isochrone with logt = 9, for which it spans AMboi = 4.4 
and logTefj = 0.22. On the observational side the corresponding values are AK = 5.5, 
A{V - K)=3.2 and A(J - K) = 1.3. 

Tables 1 to 8 list, for each age, integrated near-IR colors and their respective 
fiuctuations for three IMF slopes after the AGB to have been extended up to the end of 
the thermal pulses phase. SSPs with three values of the expected integrated V fiux are 
presented and provide an useful tool for comparison with stellar systems of different sizes. 

Indeed, when the TP-AGB is taken into account, the near-IR integrated color models 
fit satisfactorily the observed colors of LMC clusters as shown in Figure 7b. In particular, 
the logt = 9 SSP with Fy = 1000 and IMF slope x = 1.35 (Table 7) give an excelent match 
to the locus of the reddest LMC clusters in Fig. 7b. 

In order to check the sensitivity of the integrated colors on the bolometric correction 
and temperature scale adopted, that information as provided by [b'rogel fc Whitford (19877 
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for the AGB in the Galactic bulge was used to build SSPs of log t — 9. The integrated 
colors obtained at logt = 9 with the two different sets of theoretical to observational 
transformations indicated that all colors and fluctuations are smaller when the Galactic 
bulge constraints were used. The larger the stellar system, the larger the differences. For 
an SSP with Fy = 1000 the color differences are 7% in {V - K), 15% in (J - K) and 26% 
in {H -K); the color fluctuations differ by 18% in (V - K), 38% in (J - K) and 41% in 
(H-K). 



6. Concluding Renicirks 

The present work indicated the importance of stochastic effects on the near-IR 
integrated colors of young and intermediate age SSPs. If integrated observations of galaxies 
can be represented by a mixture of integrated colors of SSPs as those analyzed here, 
then it is clear that studies of stellar populations in galaxies must take into account the 
influence of bright and scarce stars on their integrated light, which can affect significantly 
determinations of age and metallicity. Gombining independent estimates (other than by 
integrated light observations) of age and metallicity for a given galaxy with its integrated 
colors would allow one, by iising the present technique, to quantify how important are the 
stochastic effects for that particular stellar population. 

We show that the inclusion of the TP- AGB phase for intermediate mass stars in the 
stellar evolution models by means of the semi-empirical approach presented here brought 
the integrated near-IR colors of intermediate age SSPs in excelent agreement with the 
spread of observed colors of LMC clusters as far as the stochastic variations are considered. 
This result reflect the expected fact that the mean integrated colors are redder and subject 
to higher fluctuations than those generated from stellar evolution models without the 
TP- AGB phase. 

In general, the main conclusions are: 

• Near-infrared bands combined to the optical V band are useful tools to the 
investigation of the red stellar population and its stochastic fluctuation in star clusters 
and, consequently, in more complex systems like galaxies. 

• The V JHK integrated colors of SSPs, when compared to those of LMC star clusters, 
gave information on the importance of taking into account TP-AGBs like Carbon 
stars in the stellar evolution models. 

• The importance of stochastic fluctuations in the integrated VJHK hght of star 
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clusters has been confirmed in this analysis, corroborating studies in the visible 
spectral range; specifically, general properties derived from the integrated light of 
composite stellar populations must account for the infiuence of stars undergoing fast 
yet bright evolutionary stages. 
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Fig. 1. — Two simulations of a young SSP with the same parameters are shown in order to 
exemphfy the stochastic effects in different CMDs. Both simulations represent SSPs having 
10^ stars and IMF slope x = 1.35. Their corresponding integrated colors are indicated. Only 
stars more massive than 1.1 Mq are plotted. 

Fig. 2. — Same as Fig. 1 for an older SSP. 

Fig. 3. — The behavior of the near-lR integrated colors and their fluctuation according to 
the number of expected post MS stars of SSPs with logt = 7.5 (a, b, c) and logt = 9.0 (d, 
e, f) are shown. Different IMF slopes are considered. 

Fig. 4. — The relative flux contribution of MS and post MS stars in VJHK to the integrated 
light of SSPs with logt = 7.5 (left) and logt = 9.0 (right). The total number of post MS 
stars is normalized to 100. The results for different IMF slopes are presented. The MS 
spectral types are indicated in the bottom panels. 

Fig. 5. — The integrated color (J — K) as a. function of SSP age. The IMF slope is x = 1.35. 
In (a) it is shown the color evolution of a massive SSP with 10® stars and a less populous 
SSP with 10^ stars. The numbers below the fluctuation bars indicate the percentage of post 
MS stars (=Npms/(100.N)). In (b) different simulation sets are presented with the number 
of post MS fixed in 100 and 500. 

Fig. 6. — The distribution of the relative VJHK flux contribution with age for upper MS 
spectral types and post MS stars, x = 1.35 and Npms=100. 

Fig. 7. — (a) The logt — 7.5 and logt — 9.0 SSP integrated colors with different number 
of post MS stars are compared to integrated colors of Large Magellanic Cloud star clusters, 
as indicated by their SWB types; the larger error bars in both SSPs correspond to the case 
where there is 1 post MS star, while the smaller ones to the case where 100 post MS stars 
are present in the system, (b) SSP integrated colors after inclusion of TP-AGB for the same 
ages as in (a); the bluer points at both ages correspond to systems with total Fy = 10 while 
the redder ones to those with total Fy = 1000. 
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Table 1. logi = 7.5 



Fv = 10 Fv = 100 Fv = 1000 







-0.5 


a; (IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


-0.5 


a;(IMF) 
1.35 


2.0 


K 




-4.567 


-5.038 


-4.423 


-7.120 


-7.038 


-6.540 


-9.542 


-9.430 


-9.201 


ctk 




4.091 


3.741 


2.575 


2.944 


2.494 


1.880 


0.659 


0.724 


0.774 


V - 


-K 


-0.467 


-0.080 


0.032 


0.534 


0.527 


0.334 


1.920 


1.680 


1.368 


(7v- 


-K 


0.894 


0.959 


0.743 


1.705 


1.500 


1.078 


0.578 


0.594 


0.653 


J - 


-K 


-0.114 


0.015 


0.064 


0.192 


0.213 


0.174 


0.731 


0.693 


0.595 




-K 


0.235 


0.245 


0.187 


0.493 


0.435 


0.324 


0.168 


0.181 


0.228 


H ■ 


-K 


-0.055 


-0.018 


-0.005 


0.017 


0.026 


0.021 


0.153 


0.146 


0.125 


(7h- 


-K 


0.061 


0.059 


0.044 


0.120 


0.101 


0.074 


0.036 


0.039 


0.053 



Table 2. logt = 7.75 



Fv = 10 Fv = 100 Fv = 1000 



a; (IMF) x(IMF) a;(IMF) 







-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


K 
ctk 




-4.422 
3.603 


-5.387 
3.617 


-4.919 
2.786 


-7.817 
2.356 


-7.371 
2.034 


-6.912 
1.769 


-9.930 
0.729 


-9.635 
0.761 


-9.334 
0.789 


V - 

(7v- 


-K 

-K 


-0.055 
1.040 


0.004 
0.966 


0.214 
0.840 


1.521 
1.557 


1.108 
1.323 


0.845 
1.155 


2.327 
0.674 


1.874 
0.690 


1.551 
0.692 


J - 


-K 

-K 


0.000 
0.281 


0.040 
0.260 


0.107 
0.202 


0.533 
0.500 


0.420 
0.432 


0.330 
0.362 


0.909 
0.301 


0.784 
0.304 


0.677 
0.312 


H ■ 

(7h- 


-K 

-K 


-0.022 
0.069 


-0.008 
0.073 


0.007 
0.048 


0.117 
0.155 


0.088 
0.127 


0.063 
0.096 


0.236 
0.131 


0.196 
0.123 


0.166 
0.120 



Tables. logi = 8 



Fv = 10 Fv = 100 Fv = 1000 



a;(IMF) x(IMF) a;(IMF) 







-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


K 




-6.133 


-4.473 


-3.985 


-7.865 


-7.154 


-6.986 


-10.490 


-10.069 


-9.718 


(7k 




3.986 


2.263 


1.703 


1.545 


1.164 


1.229 


0.772 


0.811 


0.832 


V - 


-K 


0.625 


0.416 


0.496 


1.902 


1.488 


1.249 


2.638 


2.197 


1.886 


(7v- 


-K 


1.555 


1.077 


0.833 


1.038 


0.844 


0.815 


0.708 


0.742 


0.757 


J- 


-K 


0.196 


0.155 


0.187 


0.693 


0.569 


0.496 


1.015 


0.899 


0.787 




-K 


0.430 


0.303 


0.239 


0.329 


0.274 


0.294 


0.326 


0.337 


0.349 


H- 


-K 


0.021 


0.015 


0.023 


0.133 


0.100 


0.089 


0.267 


0.228 


0.191 


(7h- 


-K 


0.106 


0.058 


0.044 


0.123 


0.083 


0.098 


0.162 


0.156 


0.157 
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Table 4. logt = 8.25 



Fv = 10 Fv = 100 Fv = 1000 



a;(IMF) x(IMF) a;(IMF) 







-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


K 




-4.293 


-5.112 


-4.011 


-7.975 


-7.520 


-7.322 


-10.616 


-10.252 


-9.907 






1.875 


2.093 


1.316 


1.585 


1.280 


1.225 


0.894 


0.862 


0.930 


V - 


-K 


0.928 


0.913 


0.816 


1.956 


1.552 


1.416 


2.830 


2.254 


1.926 


(7v- 


-K 


1.098 


1.019 


0.801 


1.057 


0.917 


0.822 


0.818 


0.795 


0.829 


J- 


-K 


0.312 


0.311 


0.284 


0.697 


0.579 


0.535 


1.085 


0.912 


0.790 




-K 


0.341 


0.289 


0.248 


0.313 


0.290 


0.278 


0.337 


0.351 


0.374 


H- 




0.044 


0.048 


0.044 


0.138 


0.110 


0.101 


0.301 


0.239 


0.200 


(7h- 


-K 


0.061 


0.056 


0.044 


0.122 


0.104 


0.104 


0.155 


0.153 


0.160 



-20- 



TableS. logi = 8.5 



Fv = 10 Fv = 100 Fv = 1000 







-0.5 


a;(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


K 




-6.886 


-5.978 


-4.502 


-8.477 


-8.137 


-7.658 


-10.977 


-10.572 


-10.335 


(7k 




3.303 


2.694 


1.424 


1.869 


1.907 


1.593 


0.804 


0.889 


0.964 


V - 


-K 


1.545 


1.241 


1.131 


2.275 


1.852 


1.630 


3.275 


2.848 


2.560 


(7v- 


-K 


1.337 


1.131 


0.813 


1.310 


1.281 


1.085 


0.752 


0.824 


0.879 


J- 


-K 


0.496 


0.392 


0.365 


0.766 


0.631 


0.561 


1.262 


1.141 


1.028 




-K 


0.340 


0.275 


0.228 


0.410 


0.381 


0.327 


0.311 


0.353 


0.382 


H- 




0.092 


0.068 


0.061 


0.171 


0.135 


0.114 


0.386 


0.342 


0.297 


(7h- 


-K 


0.101 


0.063 


0.043 


0.163 


0.138 


0.110 


0.143 


0.153 


0.161 
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Table 6. logt = 8.75 



Fv = 10 Fv = 100 Fv = 1000 







-0.5 


a;(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


K 




-6.173 


-5.472 


-5.067 


-8.541 


-8.258 


-8.222 


-11.186 


-10.773 


-10.703 






2.369 


1.890 


1.449 


1.702 


1.703 


1.627 


0.755 


0.870 


0.841 


V - 


-K 


1.626 


1.429 


1.418 


2.416 


2.087 


1.985 


3.398 


2.932 


2.773 


(7v- 


-K 


0.946 


0.831 


0.794 


1.221 


1.207 


1.144 


0.695 


0.799 


0.776 


J- 




0.503 


0.439 


0.431 


0.793 


0.688 


0.666 


1.283 


1.147 


1.112 




-K 


0.268 


0.213 


0.199 


0.413 


0.416 


0.418 


0.292 


0.362 


0.374 


H- 


-K 


0.090 


0.076 


0.074 


0.183 


0.155 


0.151 


0.400 


0.352 


0.345 


(7h- 


-K 


0.082 


0.055 


0.045 


0.171 


0.166 


0.166 


0.138 


0.164 


0.173 
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TableT. logi = 9 



Fv = 10 Fv = 100 Fv = 1000 







-0.5 


a;(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


-0.5 


x(IMF) 
1.35 


2.0 


K 




-5.613 


-4.873 


-5.719 


-8.386 


-8.041 


-8.165 


-10.863 


-10.527 


-10.581 


(7k 




1.572 


1.104 


1.573 


1.447 


1.329 


1.311 


0.680 


0.719 


0.712 


V - 


-K 


1.895 


1.719 


1.709 


2.484 


2.198 


2.161 


3.122 


2.801 


2.668 


(7v- 


-K 


0.749 


0.680 


0.638 


1.095 


0.976 


0.990 


0.633 


0.668 


0.665 


J- 


-K 


0.559 


0.503 


0.510 


0.772 


0.695 


0.689 


1.132 


1.038 


0.998 




-K 


0.188 


0.170 


0.195 


0.371 


0.352 


0.370 


0.277 


0.317 


0.328 


H- 


-K 


0.096 


0.085 


0.089 


0.171 


0.149 


0.152 


0.334 


0.303 


0.291 


(7h- 


-K 


0.058 


0.040 


0.065 


0.151 


0.143 


0.150 


0.127 


0.145 


0.149 
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Table 8. logt = 9.25 



Fv = 10 Fv = 100 Fv = 1000 



a; (IMF) x(IMF) a;(IMF) 







-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


-0.5 


1.35 


2.0 


K 
ctk 




-4.755 
0.671 


-4.851 
0.562 


-4.684 
0.603 


-7.350 
0.224 


-7.394 
0.208 


-7.252 
0.218 


-9.816 
0.081 


-9.863 
0.074 


-9.719 
0.077 


V - 

(7v- 


-K 

-K 


2.149 
0.430 


2.038 
0.373 


1.942 
0.402 


2.332 
0.165 


2.171 
0.157 


2.108 
0.167 


2.325 
0.060 


2.165 
0.057 


2.102 
0.060 


J - 


-K 

-K 


0.634 
0.128 


0.603 
0.117 


0.570 
0.130 


0.701 
0.054 


0.657 
0.054 


0.638 
0.059 


0.701 
0.020 


0.658 
0.020 


0.639 
0.022 


H ■ 

(7h- 


-K 

-K 


0.106 
0.028 


0.101 
0.025 


0.094 
0.027 


0.121 
0.014 


0.113 
0.014 


0.110 
0.015 


0.121 
0.005 


0.114 
0.005 


0.110 
0.006 
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